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Abstract 

It is shown that de-confinement can be achieved in high multiplicity non jet pp collisions at \/s= 1.8 
TeV Fermi National Accelerator Laboratory(FNAL- E735) experiment. Previously the evidence for de¬ 
confinement was the demonstrated by the constant freeze out energy density in high multiplicity events. In 
this paper we use the same data but analyze the transverse momentum spectrum in the framework of the 
clustering of color sources. The charged particle pseudorapidities densities in the range 7.0 < {dNc/drf) <26.0 
are considered. Results are presented for both thermodynamic and transport properties. The initial temper¬ 
ature and energy density are obtained and compared with the Lattice Quantum Ghromo Dynamics(LQGD) 
simulations. The energy density {e/T'^) ^ 11.5 for {dNc/drf) ^ 25.0 is close to the value for 0-10% central 
events in Au+Au collisions at y/SMN= 200 GeV. The shear viscosity to entropy density ratio(ry/s) is ~ 0.2 
at the transition temperature. The result for the trace anomaly A is in excellent agreement with LQGD 
simulations. These results confirm our earlier observation that the de-confined state of matter was created 
in high multiplicity events in pp collisions at y^=1.8 TeV. 
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1. Introduction 

The observation of high total multiplicity, high 
transverse energy, nondet, isotropic events led Van 
Hove [l| and Bjorken Q to conclude that high en¬ 
ergy density events are produced in high energy pp 
collisions [3|. These events have a far greater cross 
section than the jet production. In these events 
the transverse energy is proportional to the num¬ 
ber of low transverse momentum particles. This 
basic correspondence can be explored over a wide 
range of the charged particle pseudorapidity den¬ 
sity {dNc/dr]) in pp collisions at center of mass en¬ 
ergy ^/s = 1.8 TeV. The analysis of charged parti¬ 
cle transverse momentum data from the E735 ex¬ 
periment exhibits flow velocity of mesons and anti- 
baryons also indicating the possible evidence of 
QGP formation Q. 

Collective hydrodynamics flow has been success¬ 
ful in explaining the AA collisions at RHIC and 
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LHC energies . However pp collisions have been 
considered different from the heavy ion collisions. 
The observation of long range rapidity correlations, 
the so called “ridge”, similar to that seen in heavy 
ion collisions in high multiplicity pp collisions at 
y/s= 7 TeV by the CMS experiment suggests the 
evidence of strong radial flow d, Q . There are sev¬ 
eral theoretical papers which support the view that 
QGP can be formed in high multiplicity pp events 


In our earlier work, published in 2002, the ev¬ 
idence of hadronic de-confinement in pp collisions 
at ^/s = 1.8 TeV was presented [ll|. Based on the 
HBT analysis a constant freeze-out energy density 
^ 1 GeV/ fuA for high multiplicity pp events was 
measured [llj. The freeze-out energy density was 
found to be independent of the dNc/dij > 6. 

The objective of this work is to further analyze 
the published E735 data on the transverse mo¬ 
mentum spectra of charged particles in the frame¬ 
work of Color String Percolation Model(CSPM). 
CSPM has been successfully applied to heavy ion 
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data for thermodynamics and transport coefficients 
[iilllQE [l^ . After a brief description of 
the color string percolation model ICSPMlfl^. fl^ 
and the E735 experiment the results are presented 
for the shear viscosity to entropy density ratio(? 7 /s) 
and the Equation of State (EOS) of the deconfined 
matter. 


2. String Interactions and Percolation 

Multiparticle production at high energies is cur¬ 
rently described in terms of color strings stretched 
between the projectile and target. Hadronizing 
these strings produce the observed hadrons. The 
strings act as color sources of particles through the 
creation of qq pairs from the sea. At low energies 
only valence quarks of nucleons form strings that 
then hadronize. The number of strings grows with 
the energy and with the number of nucleons of par¬ 
ticipating nuclei. Color strings may be viewed as 
small discs in the transverse space filled with the 
color field created by colliding partons. Particles 
are produced by the Schwinger mechanisms 0- 
With growing energy and size of the colliding nu¬ 
clei the number of strings grow and start to overlap 
to form clusters 17|, [l^. At a critical density a 


macroscopic cluster appears that marks the perco¬ 
lation phase transition. 2D percolation is a non- 
thermal second order phase transition. In CSPM 
the Schwinger barrier penetration mechanism for 
particle production and the fluctuations in the as¬ 
sociated string tension due to the strong string in¬ 
teractions make it possible to define a temperature. 
Consequently the particl e sp ectrum is ’’born” with 
a thermal distribution [20|. With an increasing 
number of strings there is a progression from iso¬ 
lated individual strings to clusters and then to a 
large cluster which suddenly spans the area. In two 
dimensional percolation theory the relevant quan¬ 
tity is the dimensionless percolation density param¬ 
eter given by 00 




NSi 


( 1 ) 


where N is the number of strings formed in the colli¬ 
sions and Si is the transverse area of a single string 
and S±^ is the transverse nuclear overlap area. The 
critical cluster which spans appears for > 1-2 
110. As f increases the fraction of covered 
by this spanning cluster increases. 

We assume that a cluster of n strings behaves 
as a single string with an energy-momentum that 


corresponds to the sum of energy-momenta of the 
individual strings and with a higher color field, cor¬ 
responding to the vectorial sum of the color field of 
each individual string 00 - One can obtain the 
multiplicity /i and the mean transverse momentum 
squared ) of the particles produced by a cluster 
of n strings 0| 


l-^n — 


nSn 


-Mo; 


{pI) = 


nSi 


(P?)i 


( 2 ) 


where po and {pl)i are the mean multiplicity and 
(P() of particles produced from a single string with 
a transverse area Si = irr'i.ln the limit of high 
string density, one obtains 17. [l8| 
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where F{^) is the color suppression factor. F{^) 
is related to the 



( 4 ) 


The net effect due to F{^) is the reduction in hadron 
multiplicity and increase in the average transverse 
momentum of particles. The CSPM model calcula¬ 
tion for hadron multiplicities and momentum spec¬ 
tra were found to be in excellent agreement with 
experiment 


1^. 


It is worth noting that CSPM is a satura¬ 
tion model similar to the Color Glass Condensate 
(CGC), where (pi)i/F{0 plays the same role as the 
saturation momentum scale in the CGC model 


3. E735 experiment 


Experiment E735 was run during the 1988-1989 
Tevatron running period, primarily collecting data 
triggered to enrich high multiplicity events in p col¬ 
lisions at ^/s = 1.8 TeV 0. The E735 was located 
at the C(j, interaction region of the Fermi National 
Accelerator Laboratory (FNAL) 0, 0 271. The 
pp interaction region was surrounded by a cylin¬ 
drical drift chamber which in turn was covered by 
a single layer hodoscope including endcaps. This 
system measured the total charged particle multi¬ 
plicity 10 < A^c < 200 in the pseudorapidity range 
I 77 I <3.25. A sidearm magnetic spectrometer with 
tracking chambers and time of flight counters, pro¬ 
vided particle identified momenta spectra in the 
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range 0.1 < pt < 1.5 GeV/c. The spectrometer 
covered 0.37 < 77 < +1.00 with A4> ^ 20° {(j) is the 
azimuthal angle around the beam direction). 

The multiplicity dependence of the transverse 
momentum pt spectra was measured in pp collisions 
at y/s = 1.8 TeVji^. The invariant pt spectra was 
fitted with a power law A/{po + pt)'^ [^. Table [T] 
shows the power law ht parameters for high multi¬ 
plicity events 


27|. 


Table 1: Number of tracks Nc as measured by the E735 ex¬ 
periment in the pseudorapidity range I 77 I <3.25, {dNc/dr}), 
and the fit p arameters po and n to the invariant pt distribu¬ 
tion 


Nc 

(dNc/dp) 

Po 

n 

47 (minbias) 

7.3 

1.25 

8.35 

85 

13.07 

1.052 

7.038 

105 

16.15 

1.001 

6.743 

135 

20.76 

1.001 

6.581 

165 

25.38 

1.061 

6.766 


4 . Determination of the color suppression 
factor F(^) 


The suppression factor is determined by compar¬ 
ing the charged particle spectra from low energy pp 
collisions and high multiplicity transverse momen¬ 
tum spectra from pp. To evaluate the initial value 
of ^ from data a parameterization of pp events at 
y/s = 200 GeV is also used to compute the pt dis¬ 
tribution 


dNc/dpl = a/ipo+PtY 


( 5 ) 


where a is the normalization factor, po = 1-98 and 
n = 12.88 are parameters used to fit the data. This 
parameterization is used for high multiplicity events 
inpp collisions at y/s = 1.8 TeV to take into account 
the interactions of the strings 181. 


dNc/dp^ = 


{Poy/F{^pp)/FiYpp ) TPi) 


( 6 ) 


In pp collisions F{^pp) ~ 1 at yfs = 200 GeV due 
to the low overlap probability [I^. Figure o shows 
a plot of F{Y) as a function of charged particle mul¬ 
tiplicity per unit transverse area for high 

multiplicity events along with the minbias events. 
For pp collisions the overlap area S±_ is taken as 
the inelastic pp cross section app [iilsO,!!!!. For 


pp collision at ^ 5 = 1.8 app ^ 60 mb [^. The er¬ 
ror on F(^) is ^ 6 %. The transverse overlap area 
measured in E735 by Hanbury-Brown-Twiss is in 
excellent agreement with this inelastic cross section 


llj . The results from Au+Au collisions at y/s 


= 200 GeV [13) 01 are also shown in Fig. ([T|). 


NN 



Figure 1: The color suppression factor F{^) as a function 
of /S±{fm~^). The solid green squares are from pp at 
-v/s=1.8 TeV. The solid red circles are for Au-|-Au collisions 
at = 200 GeV(STAR data) El- The error is smaller 

than the size of the symbol. Lines are fit to pp and the STAR 
data. 


5. Temperature 

The connection between the measured F{/) and 
the temperature T{/) involves the Schwinger mech¬ 
anism (SM) for particle production. The Schwinger 
distribution for massless particles is expressed in 
terms of pj 

dn/dp^t - (7) 

where the average value of the string tension is (a:^). 
The tension of the macroscopic cluster fluctuates 
around its mean value because the chromo-electric 
field is not constant. The origin of the color string 
fluctuation is related to the stochastic picture of 
the QCD vacuum. Since the average value of the 
color field strength must vanish, it can not be con¬ 
stant but changes randomly from point to point 
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Table 2: {dNc/dr]), the measured percolation density param¬ 
eter 5, initial temperature T, initial energy density e and rj/s 
for pp -^=1.8 TeV 


{dNc/dr}) 


T (MeV) 

e{GeV/ fm^) 

T}/s 

7.30(minbias) 

0.778 

160.31 

0.50 

0.30 

13.07 

1.39 

170.73 

0.86 

0.23 

16.15 

1.42 

171.22 

1.07 

0.23 

20.76 

1.84 

178.06 

1.39 

0.21 

25.38 

2.30 

185.07 

1.75 

0.21 


[20|. Such fluctuations lead to a Gaussian distri¬ 
bution of the string tension for the cluster, which 
gives rise to the thermal distribution 12. 20l ^ 


dn / dp: 


i-pt. 


( 8 ) 


with = 'k{pI)i/F{^). The temperature is ex¬ 
pressed as 


T{0 = 



( 9 ) 


where {p1)i is the average transverse momentum 
squared of particles produced from a single string. 
At the percolation transition ^ = 1.2 and the criti¬ 
cal value of the temperature is taken = 167 MeV 
from the statistical model, which is the universal 
chemical freeze-out temperature and is a good mea¬ 
sure of the phase transition temperature [sj, 1^ . 
This determines \f{pf)i = 207 MeV [I^ . The tem¬ 
perature obtained using Eq. © for {dNc/dr]) ^ 
25 is ^ 188 MeV. This temperature is closed to ^ 
193.6 MeV obtained for 0-10 %Au+Au collisions at 
= 200 GeV[i2. 


6. Energy Density 


Among the most important and fundamental 
problems in finite-temperature QGD are the cal¬ 
culation of the bulk properties of hot QGD mat¬ 
ter and characterization of the nature of the QGD 
phase transition. The QGP according to GSPM is 
born in local thermal equilibrium because the tem¬ 
perature is determined at the string level. After the 
initial temperature T > Tc the GSPM perfect fluid 
may expand according to Bjorken boost invariant 
ID hydrodynamics |36| 




2 S±Tp 

ro 


( 10 ) 


where s is the energy density, S± the transverse 
overlap area of the colliding nuclei, and Tp^o the 
proper time. Above the critical temperature only 
massless particles are present in GSPM. To evalu¬ 
ate £ we use the charged pion multiplicity dNc/dy. 
The factor 3/2 in Eq.(12) accounts for the neu¬ 
tral pions. The average transverse mass (wt) is 
given by (to*) = \/(pt)^ + ttiq, where (pt) is the 
average transverse momentum of pion and toq the 
mass of pion. In Schwinger model Tpro is given by 

_ 2.405fi ra 

“ (mt) 

The pion multiplicity was obtained from the total 
multiplicity after subtracting the contributions due 
to kaons and protons 2^, 33- The average pt {pt) 
for pions are flat in the region 7 < {dNc/dr}) < 26. 
Table [H shows the values of e and temperature for 
various values of {dNc/drj). With the determina- 



Figure 2: e/T^ as a function of TfTc- Solid green squares 
are from pp collisions at =1.8 TeV. The lattice QCD cal¬ 
culation is shown as dotted blue line 1^ . CSPM for Au+Au 
at VTjvTv = 200 GeV are red solid circles El- 

tion of T and e one can compare the energy density 
expressed as e/T'^ with the available lattice QGD 
results Eigure[2]shows a plot of e/T"^ as a func¬ 
tion of T/Tc. The lattice QCD results are from the 
HotQCD Collaboration [38|. The result for most 
central collisions for Au-f An at ^ysNN = 200 GeV 
is also shown in Eigl2] It is observed that for high¬ 
est multiplicity {dNc/dr/) ^ 27, e/T'^ is close to the 
value obtained for Au-|-Au at 200 GeV. The mini¬ 
mum bias result is below the percolation threshold. 
This result confirms the formation of the QGP in 
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high multiplicity events in pp interactions at y/s = 
1.8 TeV. 


7. Shear viscosity to entropy density ratio 
rj/s 


In our earlier work the shear viscosity to entropy 
density ratio rj/s was obtained in the framework of 
kinetic theory and the string percolation Il4ll. The 
following expression was obtained for n/s l.lfll. 1^ 

Q. 


rj _ TL 

s 5(1 — e~^) 


( 11 ) 


where T is the temperature and L is the longitudi¬ 
nal extension of the source ^ 1 fm. 


Fig.[3]shows r]/s as a function of the temperature 
M- The lower bound shown in Fig. [3] is given by 
the AdS/CFT conjecture ^|. 

The results from Au-fAu at ^/snn = 200 GeV 
collisions are also shown in Fig. [3] for comparison 
purposes. It is seen that for {dNc/dr]) ^ 25 the p/s 
is equal to the value from the most central collisions 
in Au-f Au at 200 GeV. 



Figure 3: rj/s as a function of temperature T using Eq. mill 
[Till . Solid green squares are from pp collisions at y/s = 1.8 
TeV. CSPM results for Au-|-Au at y/Sffiy = 200 GeV are 
shown as solid red circles El- The meson gas value for rj/s 
0.7 is shown as solid black circle at T ~ 150 MeV 1^ . 
The lower bound shown is given by the AdS/CFT 0- 


8. rj/s and Trace anomaly A 


The trace anomaly (A) is the expectation value of 
the trace of the energy-momentum tensor, (0(() = 
{s — 3p), which measures the deviation from con¬ 
formal behavior and thus identifies the interaction 
still present in the medium [i^. We find that the 
reciprocal of p/s is in quantitative agreement with 
(e — 3p)/T* from LQCD over a wide range of tem¬ 
peratures [ill. Fig. 0] shows the A for pp along 
with the CSPM calculation and LQCD simulations. 
The minimum in p/s ^ 0.20 determines the peak 
of the interaction measure ~ 5 in agreement with 
the recent HotQCD values Figure 0] also shows 
the results from the Wuppertal Collaboration (45| . 
The maximum in A corresponds to the minimum in 
p/s. This is true for the highest multiplicity events 
in pp collisions. 



Figure 4: The trace anomaly A = (e—3p)/T^ vs temperature 
rmi . Solid green squares are from pp collisions at y/s = 1.8 
TeV. Blue open squares are from HotQCD Collaboration 
H . Black stars are from Wuppertal Collaboration l4fill . 


9. Equation of State: Sound velocity 

We use CSPM coupled to a ID Bjorken expan¬ 
sion. The input parameters are the initial temper¬ 
ature T, the initial energy density e, and the trace 
anomaly A are determined by data. The Bjorken 
ID expansion far the sound velocity can be written 


de " 


= ct 


5 


( 12 ) 













In the above equation the entropy density s is ex¬ 
pressed as s = {e + P)/T. The pressure P is related 
to trace anomaly P = {e — AT‘^)/3. We can express 
Cg in terms of ^ 

Cj = (-0.33) 

+ 0.0m(A/3)(^^i:^-^) (13) 
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Figure 5: The speed of sound from CSPM (red circles) and 
Lattice QCD-p4 (blue dash line) versus T/Tc[^. The green 
solid squares are from pp at =1.8 TeV. 

Since there is no direct way to obtain pressure in 
CSPM, we have made the assumption that the di¬ 
mensionless ratio A = (e — 3P)/T^ ~ l/(ry/s). 
Fig. [S] shows a plot of as a function of T/Tc. 
It is observed that the CSPM results are in very 
good agreement with the lattice calculations [38|. 
This suggests that the A as obtained from the data 
can be represented by l/{r]/s). 

10. Summary 

The pp collisions at y/s =1.8 TeV from the E735 
experiment have been analyzed in using the clus¬ 
tering of color sources phenomenology. The results 
for temperature, energy density and transport co¬ 
efficient have been extracted for high multiplicity 
events. The trace anomaly is also obtained and 
compared with the LQCD results. The clustering 


of color sources has shown that the determination 
of r]/s as a function of temperature is an impor¬ 
tant quantity that relates to the trace anomaly A. 
The main assumption of the present approach is 
that the inverse of 77 /s represents the trace anomaly, 
A = {e — 3p)/ T'*. The clustering of color sources 
(percolation) provides us with a microscopic par- 
tonic picture that also connects the transport prop¬ 
erties of the QGP to its thermodynamics. 

These results strongly suggest that even in small 
systems at high energy and high multiplicity events 
QGP formation is possible as seen in pp collisions 
at y/s = 1.8 TeV. 

The analysis of the pp collisions at y/s = 7 and 
14 TeV at the GERN LHG can map event higher 
temperatures and energy densities. Then it will be 
possible to obtain the trace anomaly at higher tem¬ 
perature as shown in Fig.(|Tl). The attenuation of 
high pt particles in high multiplicity events would 
be a second signal that the de-confined matter is 
indeed the QGP. 
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